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Abstract 

We revisit the work of Dhar and Majumdar [Phys. Rev. E 59, 6413 (1999)] 
on the hmiting distribution of the temporal mean Mf = t^^ Jq du signy^j, for 
a Gaussian Markovian process yt depending on a parameter a, which can be 
interpreted as Brownian motion in the scale of time t' = t^". This quantity, for 
short the mean 'magnetization', is simply related to the occupation time of the 
process, that is the length of time spent on one side of the origin up to time t. 
Using the fact that the intervals between sign changes of the process form a 
renewal process in the time scale t' , we determine recursively the moments of 
the mean magnetization. We also find an integral equation for the distribution 
of Mt- This allows a local analysis of this distribution in the persistence region 
{Mf — > ±1), as well as its asymptotic analysis in the regime where a is large. 
We finally put the results thus found in perspective with those obtained by 
Dhar and Majumdar by another method, based on a formalism due to Kac. 
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1 Introduction 



Consider the stochastic process yt defined by the Langevin equation 

^ = V2^e-'/\, (1.1) 

where a is a positive parameter, and rjt is a Gaussian white noise such that {rji) = 
and {Tjt^Tjt^) = 6{t2 — ti). In the new time scale 

t' = 

this process satisfies the usual Langevin equation for one- dimensional Brownian mo- 
tion, 

dyt' 



dt' 



where Ci' is still a Gaussian white noise, with {Q) = and yCf'jCt^/ = '^(^2 ~ ^i)- The 
process defined by ( |1 . 1| ) is a simple example of subordinated Brownian motion [0. 
As Brownian motion itself, it is Gaussian, Markovian, and non-stationary. 

This process appears in various situations of physical interest. For instance, it is 
described in [^] as a Markovian approximation to fractional Brownian motion. It also 
appears in ref. for the special case a = |, as describing the time evolution of the 
total magnetization of a Glauber chain undergoing phase ordering. 

Dhar and Majumdar raised the question of computing the distribution of the 
occupation time of this process, that is the length of time spent by the process on 
one side of the origin up to time t, 

where at = signyt, or equivalently of 

Mt = ] f'duau, (1.3) 
t Jo 

where Mt, the temporal mean of at, is hereafter referred to as the 'mean magnetiza- 
tion' by analogy with physical situations where at represents a spin. The distribution 
of the occupation time bears information on the statistics of persistent events of the 
process beyond that contained in the persistence exponent P, H, |^, ||, ^ |1^, |lT| . This 



exponent governs the decay ~ of the survival probability of the process, that is 
the probability that the process did not cross the origin up to time t. Actually, for 
the present case, the determination of 6 is trite, as shown by a simple reasoning [Q: 
the probability for the Brownian process yt' not to change sign up to time t' is known 
to decay as (t')~^/^, hence for the original process as This shows that 9 = a. 
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When a = |, the distribution of the fraction of time spent on one side of the 
origin by a random walker, or by Brownian motion, is given, in the long-time regime, 
by the arcsine law [0, |l|]. In contrast, when a 7^ i, the explicit determination of this 
distribution, or equivalently of the distribution of Mf, seems very difficult. However, 
as shown in Q], in the long-time regime, the computation of the asymptotic moments 
) can be done recursively, using two different methods, yielding the same results. 



The first method relies on a formalism due to Kac [|T3[ , while the second one originates 
from ref. [^. 

The method used in ref. can be applied to any (smooth) process for which 
the intervals of time between sign changes are independent, when taken on a loga- 
rithmic scale, with finite (i.e., non zero) mean i. It eventually leads to a recursive 
determination of the moments of Mf, as t — > 00 (see equation ( |3.9| ) below). 

Dhar and Majumdar make the observation that, since relations ( |3.9| ) are indepen- 
dent of i, they can be applied to the determination of the moments of Mt for the 
process (|1 . 1| ) . Comparing the resulting expressions of the moments thus obtained to 
those derived by their alternative method shows that this is indeed the case. 

However it is not obvious to understand why relations ( |3.9| ) hold for the (non- 
smooth) process ( |1 . 1| ) , since the assumptions made in order to derive them do not 
hold for such a process. In particular while, for the class of models with independent 
time intervals on a logarithmic scale, and finite i (for which the method of ref. 
has been devised), it is natural to work in a logarithmic scale of time, since the mean 
number of sign changes between and t scales as (Nf) ^ {\nt)/i, this is not so in the 
present case, since i asymptotically vanishes, and the mean number of sign changes 
scales as (Nt) ~ 2 7r~^/^t" The validity of relations ( p.9| ) for the process ( |1.1| ) 

therefore requires an explanation. 

In the present work, we revisit and extend the study done in [Q. 

We first give a new derivation of the asymptotic expressions of the moments 
I^M^^. We start from the same premises as in ref. [^, then follow another route 
-more adapted to the process under study- because of the difficulties encountered in 
applying step by step the method of ref. P to the present case (sections 2-5). We 
then identify the symmetry properties of the distributions of the random variables that 
appear in the computations, and derive a functional integral equation, the solution 
of which yields the distribution of Mt (section 6). This approach is first checked on 
the case a = \ (section 7). It is then successively applied to the study of the local 
behavior of this distribution in the persistence region, for general a (section 8), and 
to the large-a regime (section 9). 

We finally discuss some aspects of ref. We explain why a formal application 
of the method of ref. ||^ to the present case is only heuristic, and give a new inter- 
pretation of the results obtained in with the method of Kac, under the light of the 
present work (section 10). 
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2 Observables of interest 



Changes of sign of the process yt (or zero crossings) occur at discrete instants of time 
ti,t2, ■ ■ ■ ,tn, ■ ■ ■, once the process is suitably regularized at short times. We assume 
that the process starts at the origin, so that to = is also a sign change. Let Nt be 
the number of sign changes which occurred between and t, i.e., Nt is the random 
variable for the largest n for which tn<t. 

In the scale t' , where the process is (regularized) Brownian motion, sign changes 
occur at the instants of time[] t'^ = (tn)^", and Nf = Nt is the random variable 
for the largest n for which t'^ < t' . The intervals of time between sign changes are 
denoted by = t'^ — t'^_i. These are independent, identically distributed random 
variables, with a probability density function p(r'). For large values of r', p(t') decays 
proportionally to (r')^^/^. This behavior is independent of the regularizing procedure, 
while its prefactor just reflects the choice of time units. The density p(t') is therefore 
in the basin of attraction of a Levy law of index |. We choose units so that we have 
in Laplace space 



The process formed by the independent intervals of time t[,T2,---, is known as a 
renewal process. In the original scale t, the intervals of time Tn = tn — tn-i are not 
independent. 

We denote by ^at the instant of the last change of sign of the process before time t. 
This random variable depends implicitly on time t through Nf. In the scale t', we 
have 4 = (tTv)^". 

The occupation times and Tf (see equation (p!^)) are the lengths of time 
spent by the sign process at respectively in the + and — states, up to time t, hence 
t = + Tt~ . They are simply related to the mean magnetization ( p..3|) by 




(2.1) 



tMt = T+ - Tf = 2T+ ~t = t- 2T~. 



Assume that ?/t > at t = 0+, i.e., (Tt=o = +1- Then 



tMt 



( 



-it - tjv) + {tN - tN-i) if Nt = 2k + 1 (i.e., at = - 

{t - In) - {tN - tN-i) + ■■■ ifNt = 2k (i.e., at = +1). 



1), 



The converse holds if at=o = — 1- Hence we have, with equal probabilities. 



M, = ±(l-26), 



(2.2) 




Hereafter we denote by a prime any temporal variable in this scale. 
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is the fraction of time spent in the state + if at = —1, and conversely. The last 
formula can be rewritten as 

6 = Y^N, (2.3) 

where the X^r obey the recursion 

X^ = l-^X^_i, (2.4) 

with Xi = 1. Both random variables X^r and t]^_i/tis[ depend implicitly on time t 
through Nt. 

For instance, if at=Q = +1 and Nt = 4, then 

Mt = j{{t-U)-{U-h) + {t3-t2)-it2-h)+h 

= 1-26, 

where ^ = (t^ — + t2 — ti)/t = tiX^/t, with 

X, = l-^fl-^fl-^ 



3 Methods of solution 

Equations (|2.2|), ( |2.3| ) and ( |2.4| ) contain in essence the solution to the problem posed, 
namely the determination of the limiting distribution of the mean magnetization Mt 
for t ^ oo. Unfortunately, no explicit solution can be attained in general. 

However, from ( p.2|) - (|2.4|) , one can obtain recursively the moments of Mt, in the 
long-time limit. This can be done either along the lines of ref. 0, as done in 0], or by 
the method of the present work. In this section, we explain the difficulty encountered 
when applying the method of ref. to the process (|1.1[), in order to justify the more 
lengthy path we have adopted for the derivation of the moments. We shall come back 
to the comparison between the two methods in section 10. 



3.1 General framework 

Assume that, in the long-time regime, the dimensionless random variables ^Ar/t, 
t]\f_i/t]\f, Xtv, 6, and Mt possess a limiting joint distribution. Define 

H = lim — , F = lim ^ ^ , 

X = lim Xtv, ^ = lim 6, M = lim Mt. 

t— >0O t— >00 t—fOO 
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Then the equations to be solved are: 



X = I- FX, (3.1) 
e = HX, (3.2) 
M = ±(1-20- (3.3) 

These equalities hold in distribution, and the random variables entering them are not 
independent a priori. Equation ( |3.1| ) is to be understood as the fixed-point equation 
corresponding to the recursion ( p.4| ), while ( p.2| ) and ( |3.3| ) respectively correspond to 
^ and (U). 

Assume that the distribution of the random variable F is given, and that F 



is independent of X. Even so, solving ( p.l|) is difficult in general [0, [17 
However, obtaining the moments of X recursively is easier. If furthermore H and X 
are independent and the moments of H are known, then ( p.2| ) and ( |3.3| ) determine 
the moments of M. 

3.2 The diffusion equation: a reminder 

Such a situation precisely arises in the example treated in the process yt is the 
diffusing field at a fixed point of space, evolving from random initial conditions, and 
the so-called independent-interval approximation is used |T^, |T^. In the long-time 



regime, the process is stationary in the logarithmic scale of time T = Int. As a 
consequence, the autocorrelation function of the sign process, A(|AT|) = {ctO't+at) , 



only depends on the difference of logarithmic times |T^, |T9|. 

Consider the intervals of time In between successive sign changes of the process 
in the logarithmic scale of time, En = Tn — T^-i, or 

e-^- = (3.4) 

The independent-interval approximation consists in considering the intervals as 
independent, thus defining a renewal process. The distribution of the random variable 
In can then be derived, in Laplace space, from the knowledge of the correlation 
function y4(|AT|). This distribution is found to be independent of time, because 
the process is stationary in logarithmic time. Its average, (i^) = i, is some time- 
independent positive number. Explicitly, 



with 

s 



g{s) = -{l-sA{s)), (3.6) 
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where A{s) is the Laplace transform of A(T). In particular, the moments 



h = (f') = ( C-^X) = (e-''^) = Lik) (3.7) 



V t 



N 



are independent of time. Thus from ( |3.1| ) the moments of X are determined recur- 
sively in terms of the fk (see ( |A.18|) ). 



In the long-time regime, the distribution of the backward recurrence time of the 
process in the logarithmic scale, A = T — T/v, is also independent of time. This 
logarithmic recurrence time is related to the random variable H by 



-A 



e-^ = — = H. 
t 



Its distribution in Laplace space reads 



/a(^) = = (H') = .^'^'\ (3.8) 



The random variables X and H (or A) are independent. Hence (|3.2|) and ( |3.3| ) 
determine the moments of M. A remarkable fact is that the moments thus obtained, 
which are functions of the fk and of i, become independent of i when the fk are 
expressed in terms of the A{k) = Ak, using equations (|3.5|), ( p.6| ) and ( |3.7| ). Thus 



1-3^1 + 4^2 1-3^3^ ^ ^ 

^) - TT^IT ^'-'^ 

and so on. 

More generally, the method used in can be applied to any process for which the 
intervals of time between sign changes are independent, when taken on a logarithmic 
scale. It eventually leads to a recursive determination of the moments of M, resulting 
in (Kl 



3.3 The case of the process 

The situation for the process ( |1.1|) is more difficult because, in the limit t — > cx), 
tN_i/tN —>■ F = 1. Hence ( p.l|) no longer determines X, and furthermore {ijy) — > 0. 
Now, for the class of models with independent time intervals on a logarithmic 
scale, and finite (i.e., non-zero) (iiy) = ^, for which the method of ref. P|, sketched 
above, has been devised, the mean number of sign changes between and t scales as 
(Nt) ~ {\nt)/i. So, in contrast, in the present case there is no obvious reason to work 
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with this logarithmic scale of time, since (^at) asymptotically vanishes, and the mean 
number of sign changes scales as (Nt) ^ 2 7r~^/^t" [ pJ]] . 

On the other hand, if time is kept finite, then the time intervals are not 
independent and the process (|Ll]) is not stationary, precluding again the application 
of the method of ref. 

A way out is to formally apply this method to the process (|1.1|), without paying 
attention to the difficulties mentioned above, and taking advantage of the fact that 
the moments (m^^, given by ( |3.9| ), are independent on and therefore (hopefully) 
insensitive to the fact that {In) 0. This approach, which is the one followed by 
Dhar and Majumdar [Q, is however only heuristic, as further discussed in section |10[ 

Our approach relies instead on the fact that the time intervals between two 
sign changes of the process ( p. . 1| ) form a renewal process (the are independent, 
identically distributed random variables with density p(t'), given by ( p.l| ) for large 
r'). This is a fundamental property of the process (|1.1|) , and in particular of Brownian 
motion if a = |. 

This property allows us to determine the limiting distribution of H when 
t —>■ oo, and the moments fk^t of the random variable tTv-i/^Af? which are now time 
dependent. We also find the explicit time-dependent expression of (^tv)- Using the 
original equation (p.4|), instead of ( |3.1| ), and equations ( p.2|) and (|3.3|), we eventually 
recover the expressions (3^) of the^M'^^, thus extending their range of applicability. 

We then establish an integral equation for fx, and study its consequences. 



4 Distribution of 1^/1 

Using the independence of the t[,T2,---, we first determine the distribution of the 
random variable t'^, from which we then deduce that of ^at. The method used below 
is borrowed from ref. [|n|], where a thorough study of the statistics of the occupation 
time of renewal processes can be found. 

We denote by ft' the joint probability density of the random variables and A^^. 
It reads 

ft'^,Nit';y,n) = {t'^ < y , Nf = n) 

= (5(2/-4)A4<i'<d)), 

where I{t'^ < t' < t'^_^_i) = 1 if the event inside the parenthesis occurs, and if not. 
The brackets denote the average over r(, Tj, . . . Summing over n gives the distribution 
of t'^, 

oo 

ft'Jt';y) = J2ft'^Mi'-^y,ri) = {5iy-t'j,)). 

n=0 
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In Laplace space, where s is conjugate to f and u to y, 

C ft',,N{t';y,n) = 4,;v(s;«,n) = (e-< ^"^^rft'e-^*'' 
t',y \ Jt'^ I 

= (e "e " ) 



s 



p{s + uY- ^ [n>{}). (4.1) 



Note that setting m = in (|4.1| ) gives the distribution of Nt'. We finally obtain 



2^ 74 ,7V s; = - — — r . 



s(s + u) 



n=0 

In the long-time regime, i.e., for s and u simultaneously small, we get the scaling form 
which yields 

A, {t';y) ^ , ^ 

As a consequence, the random variable H' = limt/^oo t'^^^N possesses the limiting 
distribution 

h'ix) = ^ (4.2) 

TcJx{l — X) 

which is the arcsine law on [0, 1]. 

Using the equality t^/t = (t^/t')^^^", this last result yields immediately the dis- 
tribution of 



H = lim tN/t = (/f'^i/2- 

t— >oo 

which reads 



fnix) = — == = , 4.3 



This is the main result of this section. Let us define 



h{s, a) = {H^) = -B { — + -,-)= (4.4) 

,2a 



TT V2a 2'2/ V^y(— + 1 
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where B{a, b) = T{a)T{b)/T{a + b) is the beta function. For integer values of s, ( |4.4| ) 
gives the moments of fn, denoted by 

4") = h{k,a) = (h'' 



In the particular case a = |, corresponding to Brownian motion, the distribution 
of if = if' is the arcsine law ( [4.2|) , with moments 



(1/2) 1 , 1 (24-1)!! (2*:)! 



5 Determination of the moments 

In order to obtain recursion relations for the moments of the random variable X, 
we proceed in two steps. We first compute the moments of t'j^_i/t']y, from which we 
deduce those of tN-i/tN- The recursion relations for the (^X'''^ then emerge from 
( p.4|) . Equations (|3.2| ) and (|3.3|) finally determine the moments of M. 



5.1 Moments of tN^i/t 



N 



We first determine the probability density function of the joint variables t'^_i and t'^. 
In Laplace space, we have 



ut' ,-vt' / "+i j^i-st' 



p{s + u + v)''~^p{s + v)- — ^ {n>l), 

1^ ' (n = 0). 

s 

Summing over n gives 

oo 

* ra=0 

I-P(^) A I Pjs + v) \ 
s \ 1 - p{s + u + v) J ' 



(5.1) 



so that in particular ft' ^t' (s; m = 0, f = 0) = 1/s. 

The first moment of the random variable t'j^_^/t'j^ is obtained by considering 



/: ( ^ ) = / dv \-—\ C (e-"*--e^''*- 



t' \ fiw / Jo \ du I ^ t' 

= ^ + l^ln(l-p(.)) 

s s 

1 Ins 
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which leads to 



omitting the finite parts of the logarithms. 

This computation generahzes to higher-order moments, using the asymptotic 
form Q in We have 



which leads to 



k\ 



) ) - 1 - kh^^"^ (5.3) 



I' AT 

where h^^^'^'' is given by equation ( [4 .51) . In particular, since h^i^'^'^ = |, (|5.2|) is recov- 
ered. 



The result (|5.3| ) can be extended to non-integer values of k. We have thus 
(see (U)) 



t' ' 



(-AT 

,(1/2) 



as lim,fc^o ^1=1- Equation ( |5.3| ) can thus be rewritten as 



1 - kh^^''^ (4) . (5.4) 



As announced above, when t —>■ oo, the random variable t'j^_i/t'^ converges to 1, in 
distribution. 

The moments fk^t of t]^_i/t]^ are obtained from ( |5.4| ) as 

because h{k/2a, ^) = h{k, a) = h^^\ On the other hand, 
hence finally 

^ 1-A;4")|,. (5.5) 

t^OO 
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5.2 Moments of X 



From the recursion relation (E.4|), we have 



In the long-time regime, there is an asymptotic decouphng of the variables Xjv-i 
and tN-i/t^, so that it is legitimate to take Xjq — > X, while keeping the leading 
time dependence of /^^t, given by (|575|) . This procedure can be justified along the 
lines of ref. |jll|. Consider first the simple situation a = \. The difference between 
unity and tN^i/tN, which gives rise to the result (p. 51), is proportional to the interval 
tn = — tAT-i. This quantity has been shown in |]1T| to be, asymptotically for 
large t, independent of ^at, and distributed according to the a priori law p(r). A 
similar decoupling asymptotically takes place for generic a. 
Denoting the moments ^X'^'j by Xk, we obtain 

Xk = B{fk,tXk), (5.7) 

with xq = fo,t = 1, and where we have introduced the notation B for the linear 
binomial operator 



j=0 



As shown in the Appendix, ( ^.7|) implies the following recursion relations, accord- 
ing to the parity of /c, 

Xk{l + fk,t) = B{xk{l + fk,t)) (A; odd), (5.9) 

Xk{l-fk,t) = -B{xk{l - fk,t)) (A: even). (5.10) 

Using the expression ( |5.5| ) of fk^t, we obtain, in the limit t ^ oo, where £t — 0, 

Xk = B{xk) (A; odd), (5.11) 

kh^l^^Xk = -B{khi°'^Xk) (A; even). (5.12) 

These relations, which can be rewritten as 

= lT.(^)i-^y^J (^odd), (5.13) 



= -:t7Tw E . (-1)' J e^H, (5.14) 
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determine the stationary values of the Xk recursively. We thus obtain 

1 __1 3^^ 

= -7:3 + 



5.3 Moments of M 

For reasons similar to those exposed below equation ( |5.6| ), the random variables H 
and X (defined in the limit t — > oo) are independent. Thus, by ( |3.2| ) we have 



{e)=h'-^^x,, (5.15) 



which, together with equation ( |3.1| ), leads to a determination of the even moments of 
the mean magnetization M in terms of the Xk'- 

(m^) = ((1 - 2if) = B{2'' hP Xk) {k even). (5.16) 

We thus finally obtain 

'm2\ = 1-M"\ 



q !,(")/, (a) 

M^) = i+2/.r^-^^, 

M^) ^ l-54")-loi-16/.r^ + M!i(15^3"^+20.^^^^ 



(104°)+304"))4°^ 454") 4"^ ,,.7, 



h 



(a) 



and so on. 

For instance, if a = |, corresponding to Brownian motion, the successive even 
moments of M are equal to ^, |, ^, . . ., i.e.. 



\ I 22-'(j!)2 ^ 
which are the even moments of the arcsine law on [—1, 1] (see (|7.2| ) below). 
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6 An integral equation for the determination of /m 

The recursion relation ( |5.11|) expresses a symmetry property of the distribution fx'- 

fx{x) = fx{l-x) (6.1) 

(see Appendix). This is also obvious from ( p.l| ), since formally F = 1 in the present 
case. 

The recursion relation (|5.12|) , which can be rewritten as 

k(e) = -Bik(e)) (A; even), (6.2) 



expresses a symmetry property of the distribution f^, as we now show. First, it is 
easy to prove that 

Therefore ( |6.2| ) yields 

[e) = ka-0'-') (fceven). 



which is equivalent to the following symmetry property 

e/c(o = (i-o/?(i-o, (6.3) 

or 

</>(0 = 0(1-0, (6.4) 

introducing the function 

m = um- (6.5) 

On the other hand, as a consequence of ( p.2| ) and of the independence of H and 
X, the distribution is equal to the convolution of fn, given by ([4.31), and of fx- 

A«) = /'^/xW/„(i) =^r- fd.^^^, (6.6) 



hence 



In summary, two conditions determine the distribution fx{x)'- it obeys the sym- 
metry property (|6.1D, and the function (^(O, given by ( |6.7| ), obeys the symmetry 
property (|6.4| ). 

Once the probability density fx is known, is given by (|6.6| ). Finally, (|3.3|) , 
and ( |6.4D imply 

We explore the consequences of this general setup in the next three sections. 
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7 The case ol = \ 

This situation corresponds to Brownian motion. It is easy to check that the uniform 
distribution on [0, 1], 

fx{x) = 1, (7.1) 
solves the problem. Indeed, equations ( |OD and (|6.71 ) yield 



2 1-^ 



0(o = -ve(i-o, /c(o = -u , , 

which satisfy ( |6.3| ) and (|^). Finally, by (|6.8|), the limiting distribution of Mt is 
obtained: 

/M(m) = ^ (7.2) 

which is the arcsine law on [—1, 1]. 

All these results can be derived by more direct means, using the fact that in 
the present case the time intervals ri,r2, . . . between sign changes define a renewal 
process WL 



8 Local analysis in the persistence region 

The persistence region is defined by the condition M — > ±1, i.e., ^ ^ or ^ — 1. 
Considering ( p.6[ ) for ,^ ^ yields at once 

m) r,-{x~'')c-\ (8.1) 

provided the average (X~°) is convergent (see the comment below equation (|9.12|) ). 



As a consequence, using ( |6.5D and (|6.8| ), we obtain 

/M(m) ^ C{l-m^r-\ (8.2) 

with 

ol-2a 

C = ^(x-°\. (8.3) 

TT \ / 

The behavior of the distribution fx{x) as x ^ can be determined as well. 
Assuming fx{x) ^ Ax'^ (x — 0), and using ( |6.1| ), we obtain 

~ V TT r (7 + f ) 



5^1 TT JO 
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with ^ = 1 — ^,x = l — X. An identification with (|8.1|) , using again (|675|) and (|6] 
yields the values of 7 and A, hence 



fx{x) 



- 1/^4^ (^-")^°-^/^- (8-4) 



Let us compare the singular behavior ( p.2|) of /m in the persistence region with 
the beta law on [—1, 1] of same index: 

/b^^t'^(m) = C^''*" (1 - m2)"-\ (8.5) 

where 

v^r(a) ^ ^ 

A measure of the difference between the two distributions is provided by the enhance- 
ment factor 

^=^- (8-7) 

For a = |, the distribution /m(^) is the arcsine law ( |7.2D , which is a beta law. 
Equation (^ yields (X'^/^^ = 2, so that C = Cb^^t^ = l/vr, and E = 1. The 
estimate ( |8.4| ) also agrees with ( [7.1| ). 

For a 7^ |, the distribution /a/(^) is no longer a beta law, so that the enhancement 
factor E is non-trivial. 



9 Asymptotic analysis for large values of a 

For large values of a, the distributions fx{x), f^{^), and fA.i{m) are expected to share, 
at least qualitatively, some resemblance with the beta law (^.5|) . This observation 
suggests to set 

fx{x) ~ exp(-aS(x)), (9.1) 

with 

S{x) = S{l-x), (9.2) 

as a consequence of (|fci.lD. The function S{x) is expected to be regular, and positive, 
with a minimum at 5'(|) = 0, just as its counterpart 

^beta(^) = _ ln(4x(l - x)), (9.3) 
associated with the beta law ( ^.5|) . 
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With these hypotheses, 0(0; given by ( |6.7|) , can be estimated as follows. Setting 
X = ^ + £ with e < 1, we have fx{x) ^ e-"'^^^)-"^"^'^^) and a;^" -^^a ^ ^2a(^g2a£/c _ 
The change of integration variable from x to 2; = 2ae/^ yields 

m ^%,P{i)fx{i). (9.4) 



with 

P(0 = ^ / dz ^ ^ 

Setting ?/ = e~^, and returning to the variable x, we finally obtain 

r fi + 

provided the arguments of both Gamma functions are positive (this will indeed be 
the case). Note that (p.5| ) does not involve the parameter a any more. 
Omitting again pre-exponential factors, (|9.4] ) implies 



/M(m) ~ exp(-aS(i4^n. (9.6) 



In the regime of large a, the three distributions of interest are therefore given by a 
single function S{x). The problem then amounts to finding S{x), with the symmetry 
property (|9.2| ), and such that the corresponding function P(x), given by (|9.5|) , obey 



P{x) = P{l-x), (9.7) 

as a consequence of ( |6.4|) . The function S{x) is entirely determined by the above 
conditions. This property is more evident in the present regime than in the general 
case of section ||, because (|9.5|) is explicit, while (|6.7|) is an integral relationship. 

Let us first investigate the behavior of S{x) for x — |, i.e., m — >• 0, corresponding 
to the center of the distributions. Inserting the expansion 

S'(x) = C2 fx - -^^ + C4 fx - H 



in (|9.5|) , (|9.7| ), and expanding the Gamma functions accordingly, we obtain 

_ 2 _ 4 C(3) 

'''"h^' ~ + 3(h^ " (M2j^' 



and 



P(x) = i + f — 77 - 3 Vx - + 
^ ^ 2 V12(ln2)2 y V 27 
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To leading order, keeping only the quadratic term in S{x), we find that the bulk 
of the distributions is asymptotically given by narrow Gaussians for a large, namely 



fxix) ~ /^(x) ~ exp 



2a 



ln2 2 




fAii^) ~ exp 

o>i \ 2 mi 



a 



m 



The latter result is in agreement with the expressions ( ^.ITp of the moments of M, 
which behave as (M^) ^ (ln2)/a, (M*^) ^ 3(ln2)^/a^, and so on, for a 1. 

It is also worthwhile noticing that the beta law ( |8.5| ), ( |9.3| ) also becomes a narrow 
Gaussian for a large. We have S^'^^^{x) 



l\2 



and 



beta ( 
M I 



that the beta law misses a finite factor 21n2 ^ 1.3862 in the variance of the mean 
magnetization. 

The expression (|9.8| ) of the subleading amplitude C4, involving Riemann's zeta 
function, shows however that the function S{x) is altogether non-trivial. 



2.0 



1.5 - 



X 



^ 1.0 - 



0.5 





0.2 



0.4 



0.6 



0.8 



1.0 



Figure 1: Plot of the function S{x) characterizing the limiting distributions of the variables 
X and ^ and of the mean magnetization M in the large-a region, against x (full line), 
compared with the function S^^^'^{x) associated with the beta law (dashed line). 

Let us now turn to the behavior of S{x) deep in the tails of the distributions, i.e., 
for X — s> or 1, or m ^ ±1, corresponding to the persistence region. The general 
result (|8.2| ) shows that S{x) ~ —In a; has a logarithmic divergence as a; — >■ 0. As 
a consequence, the Gamma function in the numerator of expression ( |9.5| ) for P{x) 
becomes singular, as its argument goes to zero. Furthermore, in the same expression 
for P(l — x), the arguments of both Gamma functions tend to infinity. A careful 
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treatment of (|9.5|) yields the more complete expansions as x — 

S{x) = -\nx + So + 2J— + ---, P{x) = ./— + ..., (9.9) 

V vr V ^ 

while the constant 5*0 cannot be predicted by this local analysis. The square-root 
behavior of P{x) and its prefactor agree with the general results ( |8.1| ), ( |8.4| ). 



We have numerically determined the solution of ( |9.5| ), ( |9.7| ) over the whole range 
< X < |, obtaining thus accurate values of S{x). This approach yields in particular 
5*0 ~ —2.0410. Figure 1 shows a plot of the function S{x) thus obtained, compared 
with S^^'^^ix). 

As the amplitude C^^^^ of the beta law ( p.5|) remains of order unity, within expo- 
nential accuracy, the result ( |9.9| ) for S{x) implies that the amplitude C of the power 
law ( |8.2| ) in the persistence region, and the enhancement factor E defined in ( |8.7| ), 
blow up exponentially, as 



C ^ E ^ exp(Ga), (9.10) 

Q>1 Q:>1 



with 



G = lim(5^^*^(x) - Six)) = -5o - 21n2 ^ 0.6547. (9.11) 

In order to test the relevance of this large-a approach, we have evaluated numer- 
ically E for various values of the parameter a, and compared the results with the 
exponential law ( 9.10| ) predicted for large a. The computation of E can be done in 



(at least) two different ways. 

The first method consists in directly evaluating the limit 

E=\im ^ ' 



n— >oo 



The moments of the beta law ( |8.5| ) read 

beta _T (a + \)T (n + \ 



^/^lT {n + a + ^ 



while the true moments (M^") are determined from ( p.l3|) , ( ^.14|) , and ( ^.16|) , up to 



some maximal order, typically nmax = 100 to 150, beyond which numerical accuracy 
rapidly deteriorates, because the computation of (M^") involves alternating sums. 
The second method consists in combining (^.3[) and (|8.6|), getting thus 
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and in evaluating {X °) as 

:x-") = ((i-x) 



1 



^ T{n + a) 

2^ ■■ Xr. 



(9.12) 



r(a) n=0 

The behavior (p.4|) implies that the x„ decay as n~"~^/^, so that the term of order 



n in the above sum decays as n~^/^. Hence this sum is convergent, and truncating 
it at some order rimax brings a correction proportional to W^J^. The Xn are again 
determined from ( |5.13| ) and ( [5. 141 ). A linear extrapolation in n^J^ of the results of 
both schemes turns out to yield consistent results. We have for instance E 1.443 
for a = 1. 



2.0 




Figure 2: Logarithmic plot of the enhancement factor E in the persistence region, against 
a, evaluated numerically as described in the text (symbols). The straight line has the 



theoretical slope G of (|9.1lD 



Figure 2 shows our numerical results for the enhancement factor E, for values of 
a up to 3. The comparison with the exponential law ( p.lQ ) is convincing, in spite of 
the moderate values of a used. 



10 Revisiting the work of Dhar and Majumdar 



10.1 Using the method of section ^71 



Let us first show that the expressions ( [5.17[ ) for the even moments of M are the 
same as those obtained by using (|3.9| ), with the expression of A{s) appropriate to the 
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process ( |1.1| ), as done in 0]. 

The autocorrelation y4(|AT|) = {<7t<7t+at) of the sign process at in the loga- 
rithmic scale of time T = Int reads A{T) = (2/7r) arcsin(e~"'^') 0], with Laplace 
transform 



TT 



S 1 1 

2^ ^ 2' 2 



We notice that A{s) is related to h{s,a), defined in (|4.4|) , by 

h{s, a) = 1 — sA{s), 



(10.1) 



(10.2) 



where the right side is equal to 2g{s)/s by (|3.6|) . Using this identity, it is easy to 
check that the moments ( |3.9D obtained by the method of section |3.2| , with given 
by ( |10.1| ) for s = k integer, are identical to the moments ( p. IT] ) obtained by the 
method of the present work. 

It is, however, not possible to identify the intermediate results of both methods, 
as can be seen by comparing respectively equation (|5.5|) to equations (IsTTI) and (|3.5|) , 
and equation (|4.4| ) to equation (|3.8| ). This demonstrates the formal character of the 
application of the method of section |3]^ to the process (|1.1|) . (See also the discussion 
in section 3.3.) 



10.2 Comments on the results obtained using Kac's formal- 
ism 

A first comment is that the recursion relations for the coefficients Ck appearing in 
equations (14) of ref. 0] can be easily recognized to be identical to the recursion 
relations (|5.13|) , (|5.14|) for the Xk, by noting the correspondences 



a 



The second comment concerns the continuity conditions expressed in equations 
(13) of ref. 0]. Using ( p.0.3|) , these conditions yield, with the notations of the present 
work. 



ga(2X-l)\ ^ /ga(l-2X) 



These equations hold for a arbitrary, hence they are equivalent to ( p.l|) and (|6.3|) , 
respectively. 
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11 Summary and discussion 



In this paper we have revisited and extended the work of Dhar and Majumdar 
Besides providing a new recursive determination of the moments of the mean magneti- 
zation M, the present study leads to a functional integral equation for the distribution 
of the latter quantity. This framework allows a local analysis of this distribution, and 
of other relevant quantities, in the persistence region (M^ ±1), as well as a detailed 
investigation of the regime where a is large. 

The present work casts new light on the status of the expressions (|3.9|) for the 



moments of M. The method recalled in section p.2| , which leads to these relations, 
can be applied to any smooth process for which the intervals of time between sign 
changes are independent, on a logarithmic scale. For this class of processes {ijy) = I 
is finite (i.e., non-zero), and the mean number of sign changes between and t scales 
as {Nt) ^ (\nt)/l 

Relations ( |3.9| ) are also verified for the class of processes considered in this work. 
This was observed in (by comparing the expressions thus found to those obtained 
by another method, based on a formalism due to Kac), and justified by the absence 
of i in equations (|3.9D. Yet, as discussed in section 3.3, in the present case there is 
no obvious reason to work with a logarithmic scale of time, since (iiy) asymptotically 
vanishes, and the mean number of sign changes scales as {Nt) ^ 2 71-^/^^^ (See 
also the discussion in section 10.1.) 

Most of the effort of the present work was to provide a new derivation of ( |3.9|) for 
the class of processes (|1.1| ). Our approach relies on the fact that the time intervals 

between two sign changes of the process (|1.1|) form a renewal process (the 
are independent, identically distributed random variables). The derivation proceeds 



in two steps. First, relations ( p.l7| ) for the (M ) are obtained; then, using ( |10. 



equations ( p.l7| ) yield (p.9|). This extends the range of applicability of relations (|3.9| ). 
Note that for diffusion (in the independent-interval approximation) (see section 3.2), 
relations (|3l9|) hold but neither (|5.17|) nor ( |10.2|) do. 
We conclude by a few additional comments. 

In passing, let us mention another equivalent formulation of ( p.0.2| ), namely that 
the two-time autocorrelation of the sign process reads, with t < t', 

C{t,t')= f" dxfnix). (11.1) 







Another situation where ( |3.9|) , (|5.17|) , and (|10.2|) or ( |11.1|) hold is for the renewal 



processes considered in [|TT| (provided ^ < 1), which are yet another deformation of 
Brownian motion. 

Note that the first relation of ( ^.9| ), (M^) = Ai, holds whenever the two-time 
autocorrelation function is a scaling function of the ratio of the two times 0, while 
the first relation of ( ^.17|) , (M^) = 1 — (H), does not in general. For instance, for the 
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random acceleration problem, using results of ref. [20 , we find (M^) = Sv^/tt — 1^ 
0.653986 and 1 - (H) ^ 0.791335. 

This work also underlines the importance of the random variables X and H. The 
distribution of the latter is known exactly in the present case. This quantity, which 
is a natural one to consider for Brownian motion and more generally for renewal 
processes [|1T|, appears also in the context of phase ordering [plj . 



As mentioned in the introduction, the process ( |1.1| ) has been proposed as 
a Markovian approximation to fractional Brownian motion. Let us compare the 
expressions of (M^) for these two processes. For the present model we have (see (|5.17|) ) 



1 ^(- + - 

(M2) = l-^^2a 2^^ (^^2) 



2a 

while for fractional Brownian motion, with Holder index < /i < 1, we have 



(M^) = - dx arcsin ^ [ ■ (11-3) 



TT 



2x^ 



The correspondence between the two processes is made by identifying their persistence 
exponents: 6 = a = 1 — h. For 6 = 1/2, we have (M^) = 1/2 in both cases. For 
9 = 1, (|1L2D yields (M^) = 1-2/71^ 0.363380, while ( |ir3D yields (M^) = 1/3. For 
^ ^ 0, we have (M^) = 1 - cVd, with (|1L2|) yielding c = ^ 0.797885, and 



( |11.3| ) yielding c ~ 0.812233. The distributions of the mean magnetization for the 
two processes are therefore expected to be rather similar (for < < 1). 

Finally, let us comment on the changes in behaviour induced by letting the persis- 
tence exponent a vary, and compare the present process to other ones in this respect. 
The distribution of M shows a change in shape as a increases, the most probable 
value of the mean magnetization shifting from the edges to the center |Q . More pre- 
cisely, as shown in section 8, as long as a < 1, fuijn) diverges at m — > ±1, while 
for a > 1 it vanishes at these points (see equation ( ^.2| )). However for any arbitrary 
value of a the magnetization M remains distributed. 

This behaviour is actually generic, whenever the two-time autocorrelation func- 
tion of the process is asymptotically a function of the ratio of the two time variables 
In particular, this is so for diffusion. In the independent-interval approximation 
the persistence exponent 9{D) k, 0.1454a/D increases without bounds when the di- 
mension of space D is large |jl8|, [1^. As originally noted in 0, as long as 6* < 1 the 
density fMiiT^) diverges at the edges, while it vanishes there if ^ > 1. This was also 
emphasized in |^ , on the basis of scaling arguments, and recently confirmed by direct 
numerical computations 

In contrast, there are other processes for which the change in behaviour at ^ = 1 
is more radical. For fractional Brownian motion, 9 = 1 appears as a maximum 
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persistence exponent. For the renewal processes considered in ref. [jTT|, the mean 



magnetization possesses a non-trivial asymptotic distribution only if 6* < 1. 
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A Properties of the binomial operator B 

The aim of this Appendix is to prove the following property, used in section 5.2. 
Assume that the sequence Xk satisfies 

Xk = B{fkXk), (A.l) 

with Xq = fo = 1, and where B{xk) = J2j=o (j) i^^Y Xj (see (|57^ )). Then 



Xk{l + fk) = B{xk{l + fk)) (A; odd), (A.2) 

Xk{l-fk) = -B{xk{l-fk)) (A: even), (A.3) 

which are respectively equations ( |5.9|) and ( p.lOj ) in the text. 
A.l Basic properties 

In order to prove ( |A.2D and ( |A.3| ) we need the following auxiliary properties. 
First, B is its own inverse: 

B = B-\ (A.4) 



A combinatorial proof of this result can be found in ref. An alternative proof is 
obtained by noting that the action of B on exponential sequences Xk = reads 

B{y') = {l-y)K (A.5) 

This relation is invariant in the change of y to 1 — hence ( |A.4| ) follows. 
Then, we have the properties 

'^{xk) for all k, (A.6) 

?(xfc) for all k, (A.7) 

-B{xk) for all fc, (A.8) 

-B{xk) for all k. (A.9) 



Xk 


= B{xk) for /c even implies 


Xk = 


Xk 


= B{xk) for k odd implies 


Xk = 


Xk 


= —B{xk) for k even implies 


Xk = 


Xk 


= —B{xk) for k odd implies 


Xk = 


iving 


; the proofs, let us explicit the meaning 



Let us take the example of ( |A.6|) . By hypothesis, the sequence Xk satisfies the 
condition Xk = B{xk) for k even, with xq arbitrary, which is equivalent to saying that 

^ . (-1)^ Xj = {k even). (A. 10) 



j=0 



This recursion determines Xk for k odd in terms of the Xi with i = 0, . . . , k — 1 even: 
1 3 1 5 5 1 

Xl = -Xo, X3 = -X2 - -Xo, X5 = -Xi - -X2 + -Xq, 
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The property ( |A.6| ) states that Xk = B{xk) for k odd, or equivalently, 

= E C") i-iy X, [k odd), 

i=o \3 J 

which provides an infinite number of consistency relations amongst the Xk satisfy- 
ing ([AloD . 

Similarly, taking the example of ( |A.8D , by hypothesis we have x^ = —B{xk) for k 
even, with Xq = 0, which is is equivalent to 

2xk = -Y.{^'](-^y (fceven). (A.ll) 



3=0 



This recursion determines Xk for k even in terms of the x^ with £ = 1, . . . , /c — 1 odd: 

X2 = Xi, X4 = 2X3 — Xl, = 3x5 — Sxs + 3xi, 

The property (|A.8| ) states that Xk = —B{xk) for k odd, or equivalently, 

E f^l = {k odd). 



j=0 



which provides an infinite number of consistency relations amongst the Xk satisfy- 
ing ([CTD . 

We now prove the properties ( [A.6 -A.9). In order to do so, let us define, for a given 
sequence x^, the Laurent series 

oo oo 

F{z) = J2xk z-\ G{z) = B{xk) z~\ 

k=Q k=0 

We assume that these series are convergent for \z\ larger than some radius R. This 
happens e.g. if the Xk are bounded. 

The functions F{z) and G{z) are related to each other by 

Fiz) = -^Gil-z), (A.12) 
G{z) = ^F{l-z), (A.13) 



as we now show. We have 



hence, using ( [A.5|) , 
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G{z)= F{y). 



so that 

dy 
2TTiy y 

This integral is equal to the contribution of the pole at y = 1 — z, yielding (|A.13|) , 
from which ( [A. 12 ) follows. The symmetric form of the formulas ( [A. 12 ), ( A.l^j| ) is due 



to the property ( |A.4D . 

Proof of and { \K77\ ) 

The hypothesis in ( |A.6| ) implies F{z) + F{—z) = G{z) + G{—z), i.e., 

$(^) = -<l>(-^), (A.14) 



with, using (A. 13), 



^(z) = Fiz) - Giz) = Fiz) —Fil - z). 

z — 1 



Therefore {z — l)^{z) + z^{l ~ z) = 0, which can be rewritten, using (|A.14|) , as 

_ - 1) 
z z — 1 ' 

The function ^{z)/z is thus periodic, with unit period, and decaying at infinity, as 
we have ~ {xq — 2xi)/z'^ a priori. We conclude that $(2;) = identically, that 

is F{z) = G{z), implying the property ( |A.6|) . 



For the case where the Xk = {^^'^/ ^^e the moments of a random variable X, with 
density fx on [0,1], an alternative proof of ( |A.6|) is as follows. The hypothesis in 
( |A.6| ) expresses the property 



X") = [{1 - Xy) (fceven). 

As both random variables X and 1 — X are positive, this last condition is sufficient 
to imply fx{x) = — a;), hence (^X^\ = ^(1 — X)^^ for all k, which proves ( |A.6| ). 

The proof of the second property, ( |A.7|) , is very similar. The hypothesis in ( |A.7| ) 
implies 

$(z) _ _^{z - 1) 
z z — 1 

The function ^{z)/z is therefore periodic, with period two, and decaying at infinity, 
hence ^{z) = identically. 

Proof of and ( |^ ) 

The hypothesis in implies F{z) + F{-z) = -{G{z) + G{-z)), i.e., 

^(^) = -^(-2), (A.15) 
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with, using (|A.13|) , 



^{z) = F{z) + Giz) =F{z) + 



z-l 



-F(l-z). 



Therefore {z — l)'${z) — — z) = 0, which can be rewritten, using ( |A.15| ), as 

^(z) _ 1) 
z z — 1 

The function ^(z)/z is thus again periodic, with period two, and decaying at infinity, 
hence identically zero. The proof of the fourth property, ( |A.9| ), is very similar. 

A. 2 Proofs of equations (^3) 

Equation (A.l) implies the relations 



xk{i + fk) = j2(^){-iy fjx, (fcodd), 

Xk{l - fk) = J2\ ■) fj eyen), 

which determines the Xk recursively. We have thus 

1 



(A.16) 
(A. 17) 



Xi 



X2 



;i + /i)(i-/2)- 



Since the operator B is its own inverse, ( |A.1|) is equivalent to fkXk 
itself implies 

Xkil + fk) = j:(^)i-iyx, (A; odd), 

j=0 \J / 

Xkifk - 1) = XI ■ i^^y ^3 even). 



(A.18) 
B{xk), which 

(A.19) 
(A.20) 



Comparing ([A.16|) and ( |A.19| ) shows that 



EH(-iy ^.(l-/.) = (A: odd). 



hence, using the property ( [A.9| ), 

Xki^- fk) = -B{xk{l- fk)) (fceven). 
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which is equation ( |A.3|) . Similarly, comparing ( A. 17 ) and ( |A.20 ) shows that 



i=o 



or, using the property ( |A.6|) , 

+ =i3(xfc(l + /,)) (A; odd), 

which is equation (|A.2|). 
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